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We consider the universality class of the two-dimensional Tricritical Ising Model. The scaling form 
of the free-energy naturally leads to the definition of universal ratios of critical amplitudes which 
may have experimental relevance. We compute these universal ratios by a combined use of results 
coming from Perturbed Conformal Field Theory, Integrable Quantum Field Theory and numerical 
methods. 



An unifying principle in the study of critical phenom- 
ena goes under the name of universality . In the vicin- 
ity of a phase transition, when the correlation length is 
much larger than any microscopic scale, one can assign 
each system to a universality class, which is identified by 
its dimensionality D, the symmetry properties of the or- 
der parameters and the number of relevant fields. The 
first characteristic of a given universality class is the set of 
critical exponents, expressed in terms of algebraic expres- 
sions of the conformal dimensions of the relevant fields. 
Additional data of a universality class may be derived 
by the scaling properties of the free-energy alone. These 
data - called universal ratios - are pure numbers, ob- 
tained by taking particular combinations of various ther- 
modynamical amplitudes in such a way to cancel any de- 
pendence on the microscopic scales. Together with crit- 
ical exponents, universal ratios are ideal fingerprints of 
the universality classes. From an experimental point of 
view, there is by now a large literature on universal ratio 
measurements of various systems extending from binary 
fluids to magnetic systems and polymer conformations 
(for an extensive review on the subject, see ||). 

In recent years, due to the theoretical progress achieved 
in the study of two-dimensional models (at criticality 
by the methods of Conformal Field Theory (CFT) j§, 
and away from criticality by the approach of Perturbed 
Conformal Theories Qj), several universal quantities have 
been computed by different techniques for a large variety 
of bidimensional systems^such as the self-avoiding walks 
@,§, the Ising model the q-state Potts model glj, 

to name few. In this letter we will focus on the first de- 
termination of some universal ratios relative to the class 
of universality of the 2-D Tricritical Ising Model (TIM) 
for which very few universal quantities are known (see 
P,P])- Whereas the 3-D TIM describes, for instance, 
the universality class of an anti-ferromagnet with strong 
uniaxial anisotropy like FeCl2, its 2D version can describe 
the tricritical behaviour of a binary mixture of thin films 
of He 3 — He 4 |12fl or order-disorder transitions in ab- 
sorbed systems f|13|l (for a review on the theory of tricrit- 
ical points, see [|14fj. Hence there is an obvious interest in 
computing the amplest set of universal data for this uni- 
versality class and in testing the theoretical predictions 



versus their experimental determinations. 

In a continuum version of the TIM (which is, after all, 
a particular representative of this universality class), it 
is convenient to adopt a Landau-Ginzburg (LG) formu- 
lation based on a scalar field $>(x) with <E> 6 interaction. 
The LG approach permits to have a clear bookkeeping 
of the symmetry properties of each order parameter and 
to easily understand the phase diagram of the model, at 
least qualitatively. The class of universality of the TIM 
is then described by the LG euclidean action 



A = 



5 1$ + ff2 $ +33$ +34$ 



with the tricritical point identified by the bare condi- 
tions 31 = 32 = 33 = 34 = 0. Adopting a magnetic 
terminology, the statistical interpretation of the coupling 
constants is as follows: g\ plays the role of an external 
magnetic field h, 32 measures the displacement of the 
temperature from its critical value (T — T c ), 33 may be 
regarded as a staggered magnetic field h! and finally 34 
may be thought as a chemical potential [i for the vacancy 
density. Dimensional analysis shows that the upper criti- 
cal dimension of the model is D — 3, where tricritical ex- 
ponents are expected to have their classical values (apart 
logarithmic corrections). In two dimensions, although 
the mean field solution of the model cannot be trusted 
for the strong fluctuations of the order parameters, an 
exact solution at criticality is provided by CFT. In fact, 
the TIM is described by the second model of the unitary 
minimal series of CFT ||, with central charge equal to 



C 
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There are six primary fields, identified with the 



normal ordered composite LG fields ||15[, which close an 
algebra under the Operator Product Expansion (OPE). 
Only four of them are relevant (i.e. with conformal di- 
mension A < 1): cr = </?i= $ (Ai = e = <p 2 =■ < & 2 : 
(A 2 = i), a' = V3 =: $ 3 : (A 3 = ^) and i = ^ =■ $ 4 : 
(A4 = |). The fields e and t are even under the Z 2 spin- 
symmetry whereas a and a' are odd. There is another Z 2 
symmetry of the model (related to its self-duality), under 
which DeD^ 1 = —s, DtD^ 1 = t, whereas the magnetic 
order parameters are mapped onto their corresponding 
disorder parameters. Each of the above relevant fields 
can be used to move the TIM away from criticality (the 
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resulting phases of the model are discussed in |ll| ) . 

In order to derive the scaling form of the free energy 
and the set of universal ratios for the 2-D TIM, let us 
first normalise the two-point functions of the fields as 
(ipi(r)ipi(0)) ~ -^t- when r — > (in the perturbed CFT 
approach to the model, A; = 1). When the TIM is moved 
away from criticality by means of one (or several) of its 
relevant fields, with the resulting action A = Ac ft + 
J2 P 9pf d 2 xip p (x), a finite correlation length £ generally 
appears. Its scaling form may be written in four possible 
equivalent ways, according to which coupling constant is 
selected out as a prefactor 



£ = a(KigiY 



{K i9i )*. 



(1) 
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where a is some microscopic length scale, (f>ji = 1 _ A , 
and Li are universal homegeneous scaling functions of the 
ratios — Kj g V . . The terms K; are non-universal metric 
factors which depend on the unit chosen for measuring 
the external source <?j, alias on the particular realization 
of the universality class. Let f\g±, <?2, <?3, 34] be the singu- 
lar part of the free-energy (per unit volume) . According 
to which coupling constant is selected out as a prefactor, 
it can be parameterised in four possible equivalent ways 
as: 



f[gi, 92, 93, 94] = (K i; 



( Km 



(2) 



where Ti are scaling functions. For the Vacuum Expec- 
tation Value (VEV) of the fields ipj in the i th direction 
(i.e. for the off-critical theory finally obtained by gi 7^ 0, 
gk = 0, k ^ i), we have 



where, from 
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In a similar manner, 



for the generalized susceptibilities we have 
d 2 f 



r' 



where, from 



dgidgj 
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(4) 
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■ ; — KjK\K i i . These quanti- 
ties are obviously symmetric in the lower indices (r| 2 
and are respectively the usual specific heat and mag- 
netic susceptibility in the i th direction). Similarly, for 



the correlation length we have £j = a £0 9 t 



with 



£0 ~ K i _ i . From the above formulas, appropriate 
combinations can be found such that the non-universal 
metric factors Ki cancel out. Some of the 2-D universal 
ratios are: 
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In this letter we only consider the case i = 1,2, which cor- 
respond to the most important physical deformations of 
the model (the magnetic and the thermal ones), i.e. those 
which arc most accessible from an experimental point of 
view. For both the magnetic and thermal deformation 
there are no mixing among the conformal fields due to 
ultraviolet renormalization (ljjl^]. A complete analysis 
relative to all deformations of the TIM and the theoret- 
ical details of our approach will be published elsewhere 

©■ 

The e perturbation around the critical TIM is inte- 
grable and its behavior is governed by the Ej algebra 
p9| . Therefore the Bj 2 's in eq. (||) have been computed 
exactly in pc| . On the other hand, the a perturbation is 
non-integrable (numerical indications were discussed in 
p6|). In this case, the -Bji's have been numerically eval- 
uated in |22| by using the so-called Truncated Confor- 
mal Space Approach (TCSA) p^| . This method consists 
in diagonalizing the off-critical Hamiltonian on a cylin- 
der in a truncated conformal basis of the critical TIM 
such that an extimation of (tpj)i can be obtained from 
the knowledge of the eigenvectors (only the ground state 
eigenvector is needed for the VEV). All these calcula- 
tions can be easily performed by means of the numerical 
program of ref . |24| . 

In order to estimate the universal ratios, it is still nec- 
essary to calculate the r* fc 's. Their values can be ex- 
tracted in two different ways. The first method is purely 
numerical and of immediate use, since it consists in em- 
ploying the TCSA to compute numerically the deriva- 
tive -^{ l Pj)i (details will be found in j2^]). The second 
method is based on the fluctuation-dissipation theorem 
which permits to express the generalised susceptibilities 
as 



(10) 



) k = I d x(tpj{x)y k )\ 



where (• ■ -} c indicates the connected correlator. There- 
fore in this second approach we first need to evaluate the 
2— point correlation functions and then to perform the 
integration. For our calculation of the universal ratios, 
we have employed both methods, finding an agreement 
in their final outputs. Let us briefly discuss the second 
method. First of all, write the integral (10) in polar 



coordinate as T 



ok 



2ir J drr(ipj(r)(pk)i- Secondly, de- 
compose the integral over r into two integrals over the 
regions < r < R and r > R with R ~ £. When r < R, 
the correlation function (cpj(r)(pk(0))i can be efficiently 
evaluated by using a short-distance expansion [|17f 
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{<Pj(r)tpk(0))i 



(11) 



where the non-analytic dependence on the coupling con- 
stant is completly encoded into the VEV's, whereas the 
structure constants Cj fc (r) can be evaluated perturba- 
tively in g 



C l ik (r)=r 



2(A,-A J -A fc ) y^Cj^ (gi^ 2 ' 2 ^ 



(12) 



n=0 



For the TIM, C^ ' have been computed in (T(| whereas 
their first correction can be obtained by the formula 



d 



l z (tpi(oo)tpi(z)(pk(l)<Pj(0))cFT , (13) 



where the prime indicates a suitable infrared regulariza- 
tion of the integral. As shown in p5|| , an efficient way 
to compute the regularised integrals is through a Mellin 
transformation. Hence, the calculation of the above in- 
tegral ( |l3| ) on the conformal functions plus the knowl- 
edge of the various expectation values (<pi)i enables us to 
reach a quite accurate approximation of (ipj(r)(p k (0))i in 
the ultraviolet limit, i.e. for r < R. By choosing R ~ £, 
one can obtain an overlap between the ultraviolet and 
the infrared representations of the correlation functions. 
The latter is expressed by means of the spectral series 
of the correlators on the massive states | Aj.(#)) of the 
off-critical theory 



(ipj(x)(p k (0)) c 



n=l 



(14) 



where 



9n(r)=J^--~ (01^(0)14^(00 ' • • A an (0n)) x 
{A ai (0 1 )-..A a „(^)|^(O)|O) e - r SL 1 ™^- h ^ . 

As tested in several examples (see, for instance 
P,^[ [l7|j25| , |26[ | ) , the above series ( |l4| ) converges very fast 
even for r ~ £ so that its truncation to the lowest terms is 
able to capture the correct behaviour of the correlator in 
the interval r > £. For the integrable theory defined by 
the thermal deformation of the TIM, one can truncated 
the series up to the lowest 2-particle states, with the rel- 
ative matrix elements computed along the lines of the 
refs. 0,|l^,^,^8) . For the non-integrable theory defined 
by the magnetic deformation of the TIM, it is hard to go 
beyond the one-particle matrix elements and one has to 
be satisfied with the estimate of the correlators obtained 
by the one-particle contributions only: since this theory 
has two lowest masses with mass ratio m-z ~ 2mi cos ^ 
p6[ , in this case we have 

(<Pj(x)<Pk(P))i « - (fjflK ( mi r) + /? f 2 k K (m 2 r)) 



where Kq(x) is the modified Bessel function and the in- 
dices 1,2 refer to the first and second massive states. 
The one-particle matrix elements of this model fj — 
(0|(pj(0)|^4fc)can be also computed numerically by using 
the TCSA |g. 

Once an overlap of the short and large distance ex- 
pansions of the correlators in the region r ~ £ has been 
checked, a numerical integration of the correlators pro- 
vides the r* fe 's. An explicit test of the validity of the 
above method (with a corresponding estimate of its er- 
rors) is provided by the comparison of the values of T\ k 
(obtained by the numerical integration) with their ex- 
act determination extracted by the A— theorem sum rule, 
when this theorem applies 



1-A* 



-B. 



hi 



(15) 



This check shows that the uncertanties for r*. fe is at worst 
about 5%, better for the strongest relevant operators. 
Gathering all these results, a set of universal ratios for 
the TIM have been obtained. Some of them are exact, 
like (R c )\ k = ^A k , [R L )l k = §A k (k = 1,...,4). 
We have also computed those relative to the low and high 
temperature phase of the model (Table 1). An interesting 
universal ratio is provided in this case by the correlation 
length prefactors below and above the critical tem- 
perature (as extracted from the correlation function of 
the magnetic operator using its duality properties) 



£ + 



2 cos 



1.28557. 



(16) 



which can be inferred by the exact mass spectrum of the 
model and the parity properties of the excitations Jl6],[l9| . 

In summary, we have combined techniques coming 
from CFT, integrable models and numerical methods to 
obtain for the first time a set of universal quantities for 
the class of universality of the 2D Tricritical Ising Model. 
It would be interesting to have an experimental deter- 
mination of these quantities and a comparison with the 
theoretical predictions presented here. 
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Ri x = 3.54 


Kf 3 = -2.06 


R22 — 1 


R 24 = "I 


Ri 3 = 1.30 


R44 = 1 



Table 2: Universal ratios (R c )]k an ^ (Rc)^k 



{R c y 22 = 1.05 10-* 


(Rc)23 = 4 -85 10- 2 


(R c y 2i = 6.7 10-* 


{R c y 33 = 3.8 10-^ 


(R c yii = 2.0 10- 3 


(Rc)ll = -2.34 10- 2 


[R c f X3 = 1.79 10- 2 


(R c y 33 - = 3.4 10- 1 



Table 3: Universal ratio (R x )) for i,j = 1,2. 



{R x )\ = 3.897 10^2 


{R x f 2 + = 0.1111 


{R x )l = 0.116 


{R x )i~ = 0.040 


(R x y? = 


(R x y 2 - = 0.1111 



Table 4: Universal ratios i?£ and (Ra)) for i,j 
1,2-2+. 



R L e = 7.557 10~ 2 




= 1.0784 10" 1 


R'f- = 8.389 10- 1 


(R A y 2+ = 


{R A y 2 - = 3.918 IQ-' Z 


{R A )l + = 2.958 10- 1 


(R A y- = 8.260 10- 1 



Table 5: Universal ratios (Q2))k f° r *i 3i k = 1, 2 + , 2~ 



(Qa)£ +1 = 1-260 


(Qa)*--, = 1-884 


(Q 2 ) 2+2+ = 1.973 


(Q 2 )^ +2 - = 1-320 


(Q2)lt = 1-56 


{Q 2 y n = 0.442 


{Q2Y+- - 1.70 
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